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action of electromagnetic radiation is derived. The particle is neutral in charge. Equation 
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1. Introduction 

Relativistic equation of motion for perfectly absorbing spherical dust particle under the 
action of electromagnetic radiation was derived by Robertson (1937). Relativistic gener- 
alization for the case when scattered radiation is in the direction of the incident radiation 
- in proper reference frame of the particle - was presented by Klacka (1992, 2000). 

However, real particles scatter radiation in a more complicated manner. The conse- 
quence of this reality may be a completely different orbital evolution of dust particles. 
Kocifaj and Klacka (1999) show this for really shaped stationary rapidly rotating particle. 

Equation of motion for real (cosmic) dust particle was presented in Klacka (1994a) 
(and applied to real situation in Klacka 1994b). Paper by Klacka and Kocifaj (1994) 
expresses the equation in terms of optical properties used in optics. However, the last 
two papers express the result only to the first order in v/c (higher orders are neglected) 
where v is velocity of the particle, c is the speed of light. Such an accuracy may be 
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sufficient in many applications in practice. However, to be sure that the equation of 
motion presented in the above two papers is correct, one has to derive relativistically 

correct equation of motion taking into account all orders in v j c. 

The aim of this paper is to derive relativistically covariant equation of motion for real 
dust particle under the action of electromagnetic radiation. 

2. Proper reference frame of the particle — stationary particle 

The term "stationary particle" will denote particle which does not move in a given inertial 
frame of reference, although we admit its rotational motion around an axis of rotation 
(with negligible rotational velocity). Primed quantities will denote quantities measured 
in the proper reference frame of the particle. 

The flux density of photons scattered into an elementary solid angle dO' = 
sini?' dd' dip' is proportional to p'{'&' , (f') dO,', where p'{'&', if') is "phase function". Phase 
function depends on orientation of the particle with respect to the direction of the inci- 
dent radiation and on the particle characteristics; angles (p' correspond to the direction 
(and orientation) of travel of the scattered radiation, i?' is polar angle and it equals zero 
for the case of the travel of the ray in the orientation identical with the unit vector Si 
of the incident radiation. The phase function fulfills the condition 

/ p'{^',^')dn' = i. (1) 

J in 

The momentum of the incident beam of photons which is lost in the process of its 
interaction with the particle is proportional to the cross-section Cg^.^ (extinction). The 
part proportional to C'^^^ (absorption) is completely lost and the part proportional to 
(^'ext ~ ^'abs ~ ^sca (scattering) is again reemitted. 

The momentum (per unit time) of the scattered photons into an elementary solid 
angle dQ' is 

dp' sea = ls' p'{§', if') k' da' , (2) 

where unit vector in the direction of scattering is 

k' = cos-!?' Si + sini?' cos if' e\ + sini?' sinip' e-i . (3) 

S' is the flux density of radiation energy. The system of unit vectors used on the RHS 
of the last equation forms an orthogonal basis. The total momentum (per unit time) of 
the scattered photons is 

p'sca = - s' c',,, [ p'{^',^')K' dn' . (4) 
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The momentum (per unit time) obtained by the particle due to the interaction with 
radiation is 



(5) 



As for the energy, we suppose that it is conserved: the energy (per unit time) of the 
incoming radiation E-, equals to the energy (per unit time) of the outgoing radiation 
(after interaction with the particle) E'^. We will use the fact that time t' = t, where r is 
proper time. 

For the sake of brevity, wc will use "effective factors" Q'^xx instead of cfFcctivc cross- 
sections C'^xx- ^'xxx — Q'xxx ^ 1 where A! is geometrical cross-section of a sphere of 
volume equal to the volume of the particle. Equation (5) can be rewritten to the form 



dp' _ 

[ - < sint?' cos<^' > g;j ei' + [ - < sini?' sin^' > Q'^J 62'} , (6) 



^ ^ \S' A! {[O;,, - < cost?' > Q',,J Sl + 



or, in a short form 

^ = - 5' A' |g'« si + Q'l e'l' + 62'} . (7) 
ore'. J 

2.1. Summary of the important equations 

Using the text concerning energy below Eq. (5) and the last Eq. (7), wc may describe the 
total process of interaction in the form of the following equations (energies and momenta 
per unit time): 

= Ei = A' S' , 
K = {1 - Q'n) p'i - (Q'l ei' + 62') E'Jc , 

p[ = {E'Jc) Si' , (8) 

The index "i" represents the incoming radiation, beam of photons, the index "o" repre- 
sents the outgoing radiation. 

The changes of energy and momentum of the particle due to the interaction with 
electromagnetic radiation are 

dE' 

ItV 
dp' 
~dV 



E'i - E', = 0, 

Pi-P'o- (9) 



As for the condition for energy in our Eq. (8), it is equivalent to Eq. (121) in Klacka 
(1992), as for the condition for momentum p'^ in our Eq. (8), it is a generalization of Eq. 
(122) in Klacka (1992). 
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3. Stationary frame of reference 

By the term "stationary frame of reference" we mean a frame of reference in which 
particle moves with a velocity vector v = v{t). The physical quantities measured in the 
stationary frame of reference will be denoted by unprimed symbols. 

Our aim is to derive equation of motion for the particle in the stationary frame of 
reference. We will use the fact that we know this equation in the proper frame of reference 
- see Eqs. (8) and (9). We have to use Lorentz transformation for the purpose of making 
transformation from proper frame of reference to stationary frame of reference. 

If we have a four- vector = {A^, A), where A^ is its time component and A is its 
spatial component, generalized special Lorentz transformation yields 

A' = A + [{-y - l)v A/v'^ - 7 A°/c] v . (10) 
The inverse generalized special Lorentz transformation is 
AO = 7 (AO' + V ■ A'/c) , 

A = A' + [{-f - 1) V ■ A'/v'^ + 7 A"' /c] V . (11) 
The 7 factor is given by the well-known relation 



7 = 1/^1 - t;2/c2 . (12) 

As for four- vectors we immediately introduce four-momentum: 

p^ = {p°,p)^{E/c,p) . (13) 

3.1. Incident radiation 

Applying Eqs. (11) and (13) to quantity {E[/c,p^) (four-momentum per unit time - 
proper time is a scalar quantity) and taking into account also Eqs. (8), we can write 



Ei = E[y{l + V /c) , 
c 



Pi = ^ {Si + [(7 - 1) v-S-/v^ + 7/c] v) . (14) 
Using the fact that p^^ = {h v,h v Si) for photon, Lorentz transformation yields 



V = V w 
Si 



' = - {Si + f(7 - 1) v-Si/v^ - 7/cl v] , (15) 
Ui L L IS 



W 

where abbreviation 
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u; = 7 (1 - V- S^/c) (16) 
is used. 

Inserting the second of Eqs. (15) into Eqs. (14), one obtains 

= (l/w) E[ , 

p, = {l/w) {E'Jc) Si . (17) 
We have four-vector = {Ei/c,pi) = (1, Si)Ei/c = {l/w, Si/w)w Ei/c = b'^ w Ei/c. 

3.2. Outgoing radiation 

The situation is analogous to the situation in the preceding subsection. It is only a little 
more algebraically complicated, since radiation may spread out also in directions given 
by unit vectors ei, 62. The relations between ei' and e"i, £2 and e"2, are analogous to 
that presented by the second of Eqs. (15): 



, 1 



e,-- = — {ij + [{^ - l)v ej/v"" - 7/c] v] , j = 1,2 , (18) 

"3 



-'J — 

W4 



where 

Wj=^{l - v-Bj/c), j = l,2. (19) 

Using Eqs. (11) and (13) to quantity {E'g/c,p'g) (four-momentum per unit time - 
proper time is a scalar quantity), we can write 



Eo = 7(^0 + v-p',) , 

V . (20) 



Po = P'o 



E' 

(7 - 1) t; • p'Jv^ + 7 ^ 



Using also = E'^ Si /c, Eqs. (8), (18), (20) and the second of Eqs. (15) yield 

E„ = Q'j^w Eij + {1 - Q'^) Ei + 

w Ei (Qi + Q'^)^ - wEi {Q'Jwr + Q'Jw2) , 



Po = (1 - Q'H)^Si + Q'^^iv - ^ q;.^ (ce,M - 7«) • (21) 



3.3. Equation of motion 

In analogy with Eqs. (9), we have for the changes of energy and momentum of the particle 
due to the interaction with electromagnetic radiation 

dE 

—, — = Ei — Eo , 
d T 
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Putting Eqs. (21) into Eqs. (22), using also Pi = (Ei/c) Si, one easily obtains 

2 

= Q'j, {Ei - w Eij) + J2 Q'jW Ei {1/wj - 7) , 

2 



dE 



dT 
d p 



= Q'r 



E, 



Si - 



w Ei 



-yv] + ^ Q 



dr \ c J '■' c 

Equations (23) may be rewritten in terms of four- vectors: 



dp^' 
~d7 



= oi (j^ - ^ x") + E o; 



where is four-vector of the particle of mass m 



four- vector of the world-velocity of the particle is 
= (7 c, 7 v) . 

We have also found two new four- vectors 
b^j = {l/wj,ej/wj) , j = 1,2 . 
It can be easily verified that: 



(23) 



(24) 



(25) 



(26) 



(27) 



i) the quantity w Ei is a scalar quantity - see first of Eqs. (17); 

ii) Eq. (24) reduces to Eq. (7) and to the first of Eqs. (9) for the case of proper inertial 
frame of reference of the particle; 

iii) Eq. (24) yields d m/d r = 0. 



4. Conclusion 

We have derived equation of motion for real dust particle under the action of electro- 
magnetic radiation. It is supposed that equation of motion is represented by Eqs. (8) 
and (9) in the proper frame of reference of the particle. The final covariant form is rep- 
resented by Eq. (24), or, using the relations Ei = w S A' (see Eq. (38) in Klacka 1992) 
and Pi = {Ei/c) Si 

dpf" S / 



dT 



Q'^ (c 6f - u^) + J2 Q'j 



(28) 



(four- vector 6f is defined below Eq. (17)). 

Within the accuracy to the first order in v/c, Eq. (28) yields 

dv^ _ SA' 
dt m c 
2 



{q^ [(1 - v-Si/c) S, - v/c 



^ Q'j (1 - 2v Si/c + V ■ Bj/cj Bj - v/c 



(29) 
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As for practical applications, the terms vjc standing at Q'^ are negligible for majority 
of real particles. (We want to stress that values of Q'— coefficients depend on particle's 
orientation with respect to the incident radiation - their values are time dependent.) 

Acknowledgements. The author wants to thank to H. Kimura, H. Okamoto and T. 
Mukai for their remark that for the purpose of light scattering theory the transformation 
p'{-d',(p') — > Cj,^^^ dCg^^/dft' is required, where C'^^^ denotes scattering cross section and 
dCg^^/dil' is the differential scattering cross section. (June 2001; see, e. g., H. Kimura and 
I. Mann: 1998, Radiation pressure cross section for fluffy agregates, J. Quant. Spetrosc. 
Radiat. Transfer 60/3, 425- A38) 
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